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INTRODUCTION 


I  . 

Modern,  advanced,  complex  defense  technologies  often  require  an  unconven¬ 
tional  design  of  propulsive  devices  which  deviates  greatly  from  an  axisymmetr ic 
configuration.  These  unorthodox  propulsion  concepts  carry  design  flexibility 
to  a  far  more  extensive  domain  of  application  and  offer  some  special  fluid 
dynamic  advantages  over  those  with  a  circular-shaped  design.  Typical  examples 
of  three-dimensional  odd-shaped  propulsive  devices  can  be  found  in: 

•  The  scramjet  powered  hypersonic  aerospace  plane,  which  uses  part  of 
the  vehicle  under-surface  to  provide  additional  flow  expansion 

•  The  F-1S  Hornet  fighter  aircraft,  which  utilizes  movable  square¬ 
shaped  exhaust  ducts  to  obtain  more  lift  and  better  maneuverability 

•  A  new  cruise-missile  concept,  which  integrates  a  series  of  rectangu¬ 
lar  exhaust  nozzles  into  the  airframe  for  low  observability 

•  A  new  aircraft  concept,  which  features  non-circular  engines  with  a 
vectoring  nozzle  mounted  in  the  wings  for  maximum  propulsive  lift 

•  The  anti-satellite  miniature  vehicle,  which  adopts  bU  circumferen¬ 
tially  mounted  square  nozzles  for  enhanced  trajectory  control 

More  often,  however,  non-axisymmetric  propulsive  mechanizations  for  thrust 
vector  control  (TVC)  are  found  in  space  launch  vehicles;  e.g.,  PAMD-II  solid 
rocket  motor  (SRM)  with  swivel  ball  joint.  Inertial  Upper  Stage  (IUS)  SRMs 
with  fluid  bearing,  Minuteman  III  first  stage  with  swivel  ball  joint,  and 
Space  Shuttle  SRM  with  flexible  bearing.  Another  classic  example  of  this 
category  is  the  Titan  III  SRM  with  canted  nozzle. 

For  many  years,  the  Titan  III  120-in.  diameter,  five-segment  SRM  (SRM1205) 
has  been  a  very  reliable  booster  for  delivering  large  payloads  into  space  for 
various  missions.  The  motor  consists  of  a  forward  closure,  five  cylindrical 
segments,  aft  closure  with  canted  nozzle,  exit  cone,  and  exit  cone  extension. 
The  nozzle  is  canted  6°  from  bottom  dead  center  (BDC)  toward  top  dead  center 


7 


(TDC)  to  provide  a  thrust  vector  which  intersects  the  central  liquid-fueled 
core  vehicle  at  its  center  of  gravity.  An  enhanced  version  of  the  motor  is 
the  5-1/2  segment  motor  (T3^*D  SRM),  which  includes  an  additional  half- 
cylindrical  segment  near  the  forward  closure.  The  T34D  SRM  was  first  utilized 
in  1982  to  successfully  boost  the  IUS  vehicle  and  spacecraft  into  low  earth 
orbit;  it  has  since  replaced  the  SRM1205  as  zero  stage  Titan  III  during 
launching  of  large  spacecraft  into  orbit.  More  recently,  a  seven-segment 
Titan  IV  SRM  is  being  developed  to  meet  an  ever  increasing  payload  demand. 

This  is  an  improved  version  of  the  seven-segment  motor  (SRM1207)  statically 
tested  in  1969  for  the  Manned  Orbiting  Laboratory  (MOL)  program.  In  the  wake 
of  the  tragic  failure  of  the  Space  Shuttle  Challenger  on  28  January  1986,  the 
Titan  IV  SRM  has  an  even  greater  role  than  that  originally  perceived  in  the 
expendable  launch  vehicle  program  for  unmanned  flight  missions.  Figure  1 
shows  the  Titan  IV  motor  assembly  and  launch  vehicle  configuration.  Figure  2 
shows  the  Titan  IV  SRM  components.  Figure  3  illustrates  the  non-axisymmetr ic 
motor  aft  closure  geometry  and  canted  nozzle  design. 


Figure  1.  Titan  IV  Motor  Assembly  and  Launch  Vehicle  Configuration 


Figure  2.  Titan  IV  Motor  Components 


Figure  3.  Titan  IV  Motor  Aft  Closure  and  Nozzle  Configuration 
on  the  Symmetry  Plane 


9 


As  a  result  of  the  wide  application  of  asymmetric  nozzle  configurations  in 
industry,  three-dimensional  nozzle  analysis  has  received  great  attention. 
Reference  1  lists  some  recent  developments  in  this  area  for  gas-only,  on--- 
phase  flows  inside  three-dimensional  nozzles.  Gas-particle,  two-phase  flows 
inside  rocket  nozzles  were  studied  in  References  2  through  7  in  one-d  inu.-ns  iona  1 
or  two-dimensional  spaces.  Reference  8  solved  three-dimension.:!,  two-phase 
flows  inside  the  supersonic  portion  of  nozzles.  However,  a  comprehensive 
gas-particle,  two-phase,  subsonic-transonic  flow  inside  a  thre-e-d  imens  iona  1 
canted  nozzle  has  yet  to  be  addressed  in  literature.  In  fact,  for  more  than  n 
quarter  century  since  the  inception  of  the  Titan  solid  rocket  motor  program,  a 
detailed  solution  of  two-phase,  subsonic-transonic  flows  inside  the  canted 
Titan  SRM  nozzle  has  yet  to  be  obtained.  The  transonic  flow  solution  is 
recognized  in  industry  as  a  bottleneck  in  many  three-dimensional  nozzle  flow- 
studies  including  all  of  the  unorthodox  propulsive  devices  mentioned  before. 

In  this  report,  a  time-dependent  numerical  scheme  (which  has  been  used 
successfully  in  Reference  7  for  axisymmetric/two-dimens ional  solutions),  in 
conjunction  with  the  three-dimensional,  body-fitted,  coordinate  transformation 
technique  of  References  9  and  10,  is  utilized  for  the  solut:on  of  both 
gas-only,  one-phase  and  fully  coupled  gas-particle,  two-phase  subsonic- 
transonic-supersonic  flows  inside  three-dimensional  nozzles  of  arbitrary 
configurations.  Detailed  results  of  the  calculation  are  presented  for  the 
Titan  IV  canted  nozzle  flow,  even  though  the  computer  program,  which  has  been 
developed  under  this  study,  is  equally  applicable  to  other  odd-shaped, 
three-dimensional  nozzle  configurations.  The  computation  from  the  present 
time-dependent  scheme  is  carried  out  to  a  low  supersonic  nozzle  exit  plane. 
Further  expansion  of  the  flow  in  the  downstream  direction  can  be  evaluated 
from  an  efficient  steady-state,  forward  marching  scheme  presented  in 
Reference  8. 
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2.  GOVERNING  EQUATIONS 


Normalized  by  the  gas-phase  stagnation  state  corresponding  to  the  condition 
at  the  inlet  plane,  the  governing  equations  written  in  weak  conservative  form 
for  an  unsteady,  three-dimensional,  two-phase  flow  are 
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The  momentum  transfer  parameter  f ^  is  defined  as 


fi  =  VCD 
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Stokes 


is  the  particle  drag  coefficient  given  in  Reference  11  and 
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The  heat  transfer  parameter,  namely  the  particle  Nusselt  number,  is 

as  (Ref.  12) 
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and  is  defined  as  follows: 
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3.  THREE-DIMENSIONAL  TRANSFORMATIONS 


From  a  general,  three-dimensional,  arbitrary  configuration  in  the  physical 
region  (x,y,z),  the  one-to-one  differentiable  transformation  to  a  rectangular 
grid  with  a  uniform  square  mesh  in  the  computational  domain  (£,  £,  n) 
shown  in  Figure  4  can  be  accomplished  by  using  a  three-dimensional,  body- 
fitted,  coordinate  system  developed  in  References  9  and  10.  This  system 
requires  the  solution  of  three  quasi-linear ,  elliptic,  partial  differential 
equations  (PDE)  with  Dirichlet  boundary  conditions.  The  method  which  utilizes 
a  successive  over-relaxation  technique  to  solve  the  quasi-linear  elliptic  PDEs 
and  to  generate  a  body-fitted  coordinate  system  is  discussed  in  detail  in 
References  9  and  10. 


Figure  4.  Three-Dimensional  Body-Fitted  Coordinates  Transformation 
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Formally  applying  the  chain  rule  of  change  of  independent  variables  for 
Eq.  (1)  results  in  the  following  conservation  laws  in  the  (£,  n)  space: 
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and  the  Jacobian  of  transformation 


(8) 


(9) 


are  evaluated  from  numerical  differentiation  of  the  converged  solution  of  the 
body-fitted  coordinate  system.  Note  that  use  of  the  Jacobian  in  a  multi¬ 
dimensional  transformation  (Ref.  13)  eliminates  the  need  to  evaluate  second- 
order  partial  derivatives  for  the  scale  factors  of  transformation.  In  the 
physical  region  described  by  the  cylindrical  coordinate  system  (x,  r,  9), 
the  independent  variable  y  is  replaced  by  r;  z  is  replaced  by  9  in  the  scale 
factors  and  the  Jacobian  of  transformation. 
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4.  CONSIDERATION  OF  NOZZLE  CANTING 


The  geometric  transformation  mentioned  before  is  applicable  to  a  wide 
range  of  irregular,  three-dimensional  configurations  with  the  same  set  of 
coordinate  systems  for  the  motor  chamber  (hereafter  referenced  as  "motor")  and 
the  nozzle.  However,  an  added  difficulty  arises  when  the  nozzle  is  canted  or 
gimballed,  thereby  creating  a  different  reference  coordinate  system  for  the 
nozzle  from  that  for  the  motor.  In  a  canted  nozzle  flow,  the  nozzle  coordinate 
centerline  is  rotated  an  angle  (3  from  that  of  the  motor  centerline  as  shown 
in  Figure  3  for  a  Titan  IV  SRM  aft  closure  configuration.  Use  of  either  the 
motor  centerline  or  the  nozzle  centerline  as  a  fixed  reference  coordinate 
centerline  for  the  entire  flow  field  is  incorrect.  Ambiguous  interpretation 
of  flow  data  ensues  when  only  one  set  of  the  fixed  reference  coordinate 
systems  is  adopted.  In  this  regard,  it  is  obvious  that  a  one-dimensional 
analysis — based  on  an  average  area  approach — or  a  two-dimensional  study--based 
on  a  simplified,  axisymmetric  geometry — will  produce  results  of  questionable 
value  to  a  three-dimensional,  canted  nozzle  flow  investigation. 

In  this  study,  separate  coordinate  systems  are  used  for  the  motor  and 
nozzle  flow  regions.  Two  separate  three-dimensional,  body-fitted,  coordinate 
systems  utilizing  the  transformation  technique  mentioned  above  are  generated 
for  the  motor  and  the  nozzle.  The  interface  between  the  two  flow  regions  is 
taken  to  be  a  plane  passing  through  the  pivot  point  of  canting  and  perpendic¬ 
ular  to  the  nozzle  centerline  as  shown  in  Figure  3.  The  flows  in  the  two 
regions  are  solved  simultaneously.  The  treatment  of  flow  variables  at  the 
interface  needs  special  attention  and  is  discussed  later.  In  this  way, 
confusion  is  eliminated  as  to  the  use  of  the  motor  or  nozzle  centerline  as  a 
reference  coordinate  centerline  for  a  canted  nozzle  flow  field.  Precise 
interpretation  of  flow  data  can  be  obtained.  This  concept  of  multiple  grids 
for  a  canted  nozzle  flow  can  be  extended  to  a  three-dimensional,  compressible 
flow  inside  a  curved  passage. 
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5.  SOLUTION  METHOD 


The  weak  conservative  formulation  [Eq.  (8)]  is  a  hyperbolic  type  with 
respect  to  the  time  variable  in  subsonic,  transonic,  and  supersonic  flows;  it 
can  be  solved  with  the  MacCormack  finite  difference  scheme  (Ref.  14)  which  has 
been  used  successfully  for  the  asymptotic  solution  of  time-dependent,  axisym- 
metric/two-dimens ional ,  two-phase  nozzle  flows  (Ref.  7)  and  the  steady-state 
solution  of  three-dimensional,  two-phase  supersonic  nozzle  flows  (Ref.  8).  In 
the  present  study,  the  initial  condition  for  one-phase  flow  is  based  on  a  one¬ 
dimensional  isentropic  assumption,  with  the  flow  vector  set  to  a  local  incli¬ 
nation  angle  from  linear  interpolation  between  the  lower  and  upper  wall  along 
the  same  grid  line  (constant  £ ,  n  line).  The  converged,  one-phase  results 
serve  as  an  initial  guess  for  the  gas-phase  in  the  gas-particle  two-phase 
flow.  Zero  lag  (X  =  X^  =  1)  between  the  particle  and  gas  phase  is  assumed 
initially  to  start  the  two-phase  flow  calculation.  This  calculation  procedure 
and  the  inclusion  of  fourth-order  damping  terms  to  the  second-order  MacCormack 
finite  difference  scheme  are  much  the  same  as  those  in  Reference  /,  except 
with  an  added  dimension  along  the  circumferential  direction. 

The  nozzle  exit  boundary  condition  is  based  on  a  linear  extrapolation  of 
conservative  variables  [D  of  Eq.  (1)],  since  the  flow  is  assumed  to  be  super¬ 
sonic  at  the  nozzle  exit  plane,  and  any  error  generated  from  the  extrapolation 
is  not  expected  to  propagate  back  and  affect  upstream  results. 

Two  choices  of  inlet  boundary  conditions  applied  to  the  motor  entrance 
region  have  been  incorporated  into  the  computer  program.  The  first  inlet 
boundary  condition  is  obtained  from  an  axisymmetric  characteristics  formula¬ 
tion  similar  to  that  discussed  in  Reference  15  and  used  in  Reference  7.  This 
boundary  condition  applies  to  an  axisymmetric  inlet  flow  but  can  be  extended 
to  non-axisymmetr ic ,  three-dimensional  inlet  flow  (Ref.  16).  The  solution 
obtained  from  applying  this  inlet  boundary  condition  usually  converges  to  a 
fairly  smooth,  low  subsonic  inlet  condition  close  to  an  isentropic,  uniform, 
one-dimensional  inlet  flow;  it  is  suitable  for  the  inlet  flow  that  is  not 
perturbed  significantly  by  crossflow  occurring  upstream  of  the  inlet  plane. 
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The  second  inlet  boundary  condition  corresponds  to  specified  fixed  flow 
variables  at  the  inlet  plane.  For  one-phase  flow,  Mach  number,  flow  angles, 
and  static  pressure  are  specified  on  the  inlet  plane.  For  two-phase  flow, 
particle  velocity  lag  (X  )  and  gas-to-particle  temperature  ratio  (X^)  are  also 
needed  on  the  inlet  plane.  This  inlet  boundary  condition  is  suitable  for  a 
non-uniform  flow  at  the  inlet  plane,  where  flow  variables  can  be  obtained  from 
experimental  data  or  theoretical  analysis.  In  applying  this  inlet  boundary 
condition,  a  second  order  interpolation  along  the  axial  main  flow  direction  is 
desirable  to  smooth  the  flow  variables  at  the  grid  points  adjacent  to  the 
inlet  plane  in  order  to  avoid  instability  near  the  motor  inlet  region. 


For  a  cylindrical  coordinate  system,  the  conservative  variables  G  and 
in  Eq.  (1)  at  motor  and  nozzle  centerlines  are  obtained  from  a  linear  extrapo¬ 
lation  of  the  values  at  the  interior  points.  In  decoding  conservative 
variable,  D,  the  flow  variables  at  the  motor  and  nozzle  centerlines  are 
obtained  from  interpolating  the  flow  variables  at  one  grid  point  above  and 
at  one  grid  point  below  the  centerline  on  the  symmetry  plane.  In  a 
non-axi symmetric  nozzle,  the  radial  velocity  components  at  the  coordinate 
centerline  are  not  necessarily  zero.  On  the  symmetry  plane,  however,  the 
meridional  velocity  components  are  zero.  Therefore,  at  0  =  0  on  the 
symmetry  plane,  the  following  relations  are  used: 


fl,l  f 2 , 1  '  (f2,l  f2,NK)  (r 


2,1 


2,1  +  r2,NK) 


where  f  .  , 
J»k 
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stands  for  p,  u,  e,  p  ,  u ^ ,  and  h ^ ,  and  the 
indexes  along  r,  0  directions  on  a  constant 
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where  g.  ,  stands  for  v  and  v.;  and 
J ,  k  J 
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where  h.  ,  stands  for  w  and  w.. 
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The  scalar  flow  variables  and  -elocity  vectors  are  uniquely  defined  at  the 
coordinate  centerlines  from  the  equations  previously  shown.  The  components  of 
the  velocity  vector,  however,  take  different  values  at  different  meridional 
planes.  Therefore,  at  9  ^  0  on  the  symmetry  plane  centerlines,  the  flow- 
variables  are  computed  from  those  at  0  =  0  as  follows: 


g 


l,k 


g  cos  0  +  h  sin  9 

1,1  1  1  1 


h  =  -g  sin  9  +  h  cos  9 

1  ,  K  1,1  1,1 


In  this  way,  the  flow 
and  nozzle  in  cylindrical 
fashion  without  ambiguity 


variables  at  the  coordinate  centerlines  of  the  motor- 
coordinates  are  evaluated  in  a  three-dimensional 
and  undue  complication. 


The  flow  variables  at  the  wall  boundary  are  obtained  from  linear 
extrapolation  of  data  from  the  adjacent  interior  points  and  then  modified 
by  the  local  tangency  condition,  q  •  n^  =  0,  with  three  components  of  velocity 
vector  found  as  follows: 

u  =  |  q  |  nfc 

v  =  |  q  j  nt 

w  =  |  q  j  nt 


where  3  is  the  gas  velocity  vector  on  the 
in  x,  r,  0,  directions,  respectively. 


wal  1 


and 
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the  unit  vectors 
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-»  M  ^  w  w 
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is  the  unit  tangent  vector  on  the  wall  and 


is  the  unit  normal  vector  to  the  wall  for  [f(x,  r,  G ) ]  =  [r  -  R  (x,  6)]  =  0 

w  w 

in  a  cylindrical  coordinate  system  or 


for  [f(x,  y,  z )  ]  ^  =0  in  a  rectangular  coordinate  system.  Similar  expressions 
are  applied  to  the  particle-phase  velocity  components,  even  though  it  is 
recognized  that  particles  would  interact  with  the  solid  wall  in  a  complicated 
manner,  which  accounts  for  enhanced  motor  and  nozzle  wall  erosion.  The 
particle-boundary  interaction  is  a  subject  of  further  research. 


On  the  interface  plane  of  the  motor  and  nozzle  flow  regions,  the  flow 
variables  are  evaluated  from  interpolating  variables  at  upstream  and 
downstream  nodal  points  on  the  motor  and  nozzle  interface,  respectively,  as 
illustrated  in  Figure  5.  The  scalar  flow  variables  (pressure,  density, 
temperature,  and  total  energy)  on  the  interface  plane  are  invariant  with 
respect  to  coordinate  change  from  the  motor  to  nozzle  centerline  and  can  be 
obtained  directly  from  this  interpolation  process.  The  velocity  vector, 
however,  has  different  components  on  the  interface  for  the  motor  and  nozzle, 
respectively,  and  can  be  found  by  the  following  method. 


UPSTREAM 

POINTS 

INTERFACE 


DOWNSTREAM 
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MOTOR  COMPUTATIONAL  DOMAIf 


NOZZLE  COMPUTATIONAL  DOMAIN 


Figure  5.  Determination  of  Flow  Variables  at  Interfac 


Let  the  velocity  at  any  point  in  the  motor  coordinates  be 


q  =  Vxl  +  Vrl  +  W1 101 


where  i  , ,  i  i«,  are  the  unit  vectors  in  the  motor  coordinates.  This  same 
xl  rl  01 

vector  at  the  same  point  in  the  nozzle  coordinates  is 


q  =  Vx2  +  V21r2  +  V02 


where  i  i  irt„  are  the  unit  vectors  in  the  nozzle  coordinates. 
x2  r2  02 


Since  both  the  motor  and  nozzle  coordinates  are  orthogonal,  the  dot 
products 
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Note  /  a.. a,  .  =  6..  =  Kronecker  delta  function. 
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For  a  nozzle  with  canting  angle  (3,  the  direction  cosines  are 

cosB  -sinBcosG^  sinfsinG^ 

a..  =  sinBcosG,  cosBcosG, cos9„  +  sin9,sin0„  -cosBcosG, s inG„  +  sinG..cos0„ 
lj  1  1212  1212 

-sinBsinG^  -cosBs in9^cos90  +  cosG^sinQ^  cosBsin9^sin90  +  cosG^cosG,, 
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and  the  two  coordinates  are  related  by  the  following  equations: 


r0sin90 

tar"'l  x0sinB  +  r9cosBcos09 
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X2  =  x^cos^  +  r^sinBcos9^ 


(13) 


(1M 


which  also  imply  that  uniform  circumferential  grid  division  in  09  (nozzle 
coordinates)  will  result  in  non-uniform  circumferential  grid  division  in  ©^ 
(motor  coordinates)  on  the  interface  plane  for  a  nozzle  with  non-zero  canting 
angle . 


Equation  (10)  is  used  to  define  the  components  of  velocity  vectors  on  the 
nozzle  side  with  respect  to  the  motor  coordinates,  and  Equation  (11)  is  used 
to  define  the  components  of  velocity  vectors  on  the  motor  side  with  respect  to 
the  nozzle  coordinates.  The  velocity  components  at  the  interface  are  found 
from  interpolating  the  velocity  components  at  the  upstream  and  downstream  grid 
points  adjacent  to  the  interface  with  a  consistent  reference  (motor  or  nozzle) 
coordinate  system. 
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6.  THREE-DIMENSIONAL  BODY-FITTED  GRIDS 


The  initial  guess  of  three-dimensional  body-fitted  grids  for  both  motor 
and  nozzle  are  based  on  linear  interpolation  from  a  specified  boundary  nodal 
points  arrangement,  which  consists  of:  (a)  the  two-dimensional  boundary- 
fitted-coordinate  systems  (Ref.  17)  in  upper  (TDC,  0  =  0)  and  lower  (BDC,  0  =  it 
symmetry  planes  and  at  90°  (0  =  ir/2)  plane  for  both  motor  and  nozzle  flow 
regions;  (b)  a  uniform  circumferential  division  for  nozzle  and  a  non-uniform 
circumferential  division  for  motor  according  to  Eqs .  (13)  and  (1A);  and  (c)  the 
linearly  interpolated  grid  radial  coordinates  between  TDC  and  90°  and  between 
90°  and  BDC  for  motor  inlet,  interface,  and  nozzle  exit  planes.  Care  has  been 
taken  to  ensure  that  grid  points  are  continuous  on  the  coordinate  centerlines 
and  on  the  interface  plane  that  joins  the  motor  and  nozzle  flow  regions. 

These  boundary  nodal  meshes  are  presented  in  Figure  6  for  the  Titan  IV  motor 
and  nozzle  configuration  shown  in  Figure  3  when  the  aft  closure  propellant  is 
completely  consumed,  corresponding  to  a  time  approximately  78.3  sec  into  motor 
burn  for  a  nominal  motor  burn  time  of  123  sec. 


mm 

M 


•  SYM  PLANE  EXIT  PLANE 


BDC  BDC 
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Figure  6.  Boundary  Grids  on  Motor  Entrance,  Motor  Symmetry, 
Nozzle  Symmetry,  and  Nozzle  Exit  Planes 


The  inhomogeneous  terms  in  the  boundary-fitted  coordinate  transformation 
have  been  turned  on  to  keep  the  streamwise  grid  lines  near  the  nozzle  wall  in 
the  concave  region.  For  irregular  motor  and  nozzle  configurations  whose 
boundary  geometry  cannot  be  satisfactorily  described  by  linear  interpolation 
of  radial  coordinates  between  TDC  and  90°,  and  between  90°  and  BDC  planes, 
provisions  have  been  incorporated  in  the  program  to  allow  for  geometry  input 
of  boundary  radial  coordinates  on  any  specified  meridional  planes  in  addition 
to  those  at  TDC,  90°  and  BDC. 

The  three-dimensional  interior  grid  points  are  generated  from  a  successive 
over-relaxation  iteration  scheme  for  the  solution  of  the  three-dimensional 
quasi-linear  elliptic  system  mentioned  before.  With  the  careful  arrangement 
of  boundary  nodal  meshes,  fast  convergence  of  iteration  for  the  three- 
dimensional  interior  grid  generation  is  assured.  It  is  important  that  the 
inhomogeneous  terms  in  the  three-dimensional  body-fitted  coordinate  system  for 
the  interior  nodal  points  generation  be  included,  to  be  compatible  with  the 
two-dimensional  boundary-fitted  coordinates  specified  on  the  three-d imens ional 
grid  boundaries  for  nodal  points  clustering  control.  An  exponential  decay 
function,  similar  to  that  in  Reference  17,  is  applied  for  the  three-dimensional 
interior  coordinate  control.  It  takes  17.2  sec  to  generate  a  three-dimensional 
Titan  IV  motor  grid  (NI  =  AO,  NJ  =  18,  NK  =  13  in  the  x,  r,  0  directions, 
respectively)  and  1.2  sec  to  generate  a  Titan  IV  SRM  three-dimensional  nozzle 
grid  (20  x  18  x  13)  on  a  CRAY  X-MP1A  super  computer.  Numerical  differentiation 
is  then  utilized  to  calculate  the  scale  factors  of  transformation  from  an 
irregular  three-dimensional  physical  region  to  a  uniformly  divided  rectangular 
parallelepiped  computational  domain.  This  process  takes  only  0.1  sec  for  the 
motor  and  0.01  sec  for  the  nozzle  flow  region.  The  hidden-line  plots  of  the 
converged  Titan  IV  SRM  three-dimensional  grids  are  shown  in  Figure  7,  and  the 
wall  boundary  geometry  for  the  motor  and  nozzle  is  given  in  Figure  8. 
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Figure  7.  Three-Dimensional  Body-Fitted  Coordinates  Grids 
for  Motor  and  Nozzle 


Figure  8.  Motor  and  Nozzle  Wall  Geometry 
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7.  ONE-PHASE  COLD  FLOW,  y  =  1.4 


To  aid  design  and  evaluation  of  the  full-scale  Titan  IV  SRM,  a  cold  flow- 
test  with  a  4.65%  scale  model  was  conducted  (Ref.  18).  Gaseous  nitrogen 
(y  =  1.4)  was  introduced  uniformly  into  the  scaled  motor  chamber  through 
sintered  bronze  porous  surfaces  to  simulate  propellant  burning.  Nitrogen  was 
introduced  at  a  chamber  pressure  of  45  psia  and  chamber  temperature  of  492°R. 
The  model  nozzle  throat  diameter  was  1.895  in.  An  exact  scale-down  of  the 
motor  aft  closure  and  nozzle  configurations  was  utilized  in  the  cold  flow- 
test.  This  model  provided  test  data  for  comparison  with  results  from  the 
present  theoretical  calculation. 

In  the  one-phase  cold  flow  (y  =  1.4)  numerical  study,  the  fixed  non- 
uniform  inflow  condition  is  taken  from  the  three-dimensional  cold  flow  test 
data  at  TDC  and  BDC  planes  as  illustrated  in  Figure  9.  An  isentropic  one¬ 
dimensional  Mach  number  at  the  same  area  ratio  is  also  indicated  in  the  figure 
for  comparison.  Inflow  Mach  numbers  at  other  meridional  planes  are  obtained 
from  linear  interpolation  of  the  Mach  numbers  at  TDC  and  BDC  planes.  The 
computed  dimensionless  mass  flow  rate,  wall  Mach  number  at  TDC,  90°,  BDC,  and 
centerline  Mach  number  at  the  throat  station  for  every  integration  step  are 
shown  in  Figure  10.  The  initial  throat  centerline  Mach  number  is  less  than 
one,  because  the  grid  plane  (constant  £  plane)  passing  through  the  throat 
station  at  the  nozzle  wall  is  initially  taken  to  be  a  sonic  surface  which 
curves  toward  downstream  at  the  centerline  in  the  physical  region.  This  helps 
to  speed  up  the  solution  convergence.  Note  that 

W  =  f  f (pu  +  p ,u.)dA/[(2f)^p  ,V  ,L2)  (15) 

JJ  3  3  tl  maxi 

is  the  dimensionless  mass  flow  rate  evaluated  at  the  nozzle  geometric  throat. 
For  a  one-phase  flow,  the  contribution  from  particle  mass  flux  is  zero. 
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The  transonic  flow  field  in  the  vicinity  of  nozzle  throat  converges 
rapidly.  It  is  the  resolution  of  the  motor  subsonic  flow  field  which  consumes 
most  of  the  machine  calculation  time.  The  convergence  criterion  used  in  the 
transonic  flow  regime  requires  that  the  difference  in  mass  flow  rate  and  in 
Mach  number  on  TDC,  90°,  BDC,  and  centerline  at  the  nozzle  throat  be  less  than 
0.001%  for  three  consecutive  time  integration  steps.  The  vectorized  computer 
program  developed  under  this  study  with  40  x  18  x  13  grid  points  in  the  motor 
and  20  x  18  x  13  in  the  nozzle  flow  regions  takes  up  1.1  million  words  memory 
on  a  CRAY  X-MP14  super  computer.  The  converged  Titan  IV  one-phase  cold  flow 
transonic  solution  takes  2083  integration  steps  which  require  4  min,  54  sec 
computation  time.  The  computed  three-dimensional  nozzle  cold  flow  discharge 
coefficient  is  0.9855,  based  on  a  converged  three-dimensional  gas  mass  flow¬ 
rate  of  2.997  lb/sec  and  a  one-dimensional  isentropic  gas  mass  flow  rate  of 
3.041  lb/sec  at  the  nozzle  throat  for  the  4.65%  scale  cold  flow  model. 

The  Mach  number  distributions  on  the  wall  and  along  the  centerlines  are 
depicted  in  Figure  11.  It  shows  that,  as  flow  approaches  the  motor-nozzle 
interface  plane,  higher  Mach  number  occurs  on  the  TDC  wall  surface  with 
steeper  slope  change  than  that  on  the  BDC  wall  surface.  The  three- 
dimensionality  of  flow  is  evident  from  these  calculated  Mach  number  distribu¬ 
tions.  The  computed  Mach  numbers  at  the  throat  are  1.339,  1.324,  1.293,  and 
0.813  at  TDC,  90°,  BDC,  and  centerline,  respectively,  for  the  cold  flow  with 
y  =  1.4.  The  results  of  a  single  one-dimensional  analysis  at  TDC  +  90° 
plane  are  also  shown  in  the  figure  for  comparison.  It  clearly  indicates  the 
inadequacy  of  the  one-dimensional  analysis  applied  to  a  three-dimensional 
configuration.  For  the  Titan  IV  nozzle  which  is  canted  from  BDC  toward  TDC, 
the  higher  pressure  load  associated  with  lower  flow  speed  on  the  BDC  plane 
than  that  on  the  TDC  plane  of  the  nozzle  is  to  be  expected.  Shown  on  the  same 
figure  are  the  data  obtained  from  the  cold  flow  test  (Ref.  18).  In  general, 
the  results  of  the  present  three-dimensional  flow  analysis  are  in  very  good 
agreement  with  cold  flow  test  data,  except  in  the  low  speed  compression  corner 
region  of  the  aft  closure,  where  the  boundary  layer  is  thick  and  viscous 
effects  prevail. 
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Figure  11.  Gas  Mach  Number  Distribution  (one-phase,  y  =  1.4) 

Figure  12  shows  the  computed  Mach  number  contour  on  the  TDC-BDC  symmetry 
plane.  An  isometric  projection  of  the  three-dimensional  Mach  number  contour 
plot  is  given  in  Figure  13  for  a  clear  visualization  of  the  computed  overall 
three-dimensional  flow  field.  The  sonic  surface  in  the  nozzle  is  marked  in 
Figures  12  and  13.  This  sonic  condition  is  of  particular  interest  to  rocket 
nozzle  designers,  not  only  because  it  provides  a  "sonic  barrier"  beyond  which 
an  efficient  supersonic  marching  scheme  or  method  of  characteristics  is 
appli  cable  and  shock  waves  usually  emerge,  but  also  it  is  the  state  at  which 
maximum  convective  heating  occurs  at  the  nozzle  wall.  Downstream  of  the  sonic 
state,  flow  density  and  pressure  decrease  drastically,  often  resulting  in 
reduced  convective  heating  to  the  nozzle  wall. 
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Note  that  the  Mach  number  contour  curves  are  a  direct  reflection  of  the 
calculated  three-dimensional  flow  field  from  the  present  analysis  and  are 
obtained  directly  from  a  Calcomp  pen  plotter  based  on  the  information  stored 
in  magnetic  tapes.  The  smoothness  of  the  contour  curves  and  near-normal 
intersection  of  the  curves  with  centerlines  of  the  motor  and  nozzle  indicate 
that  the  numerical  scheme  used  in  the  present  analysis  for  boundary  treatment 
works  well.  Similar  to  that  discussed  in  Reference  7,  a  small  recompression 
zone  occurs  at  the  wall  in  the  supersonic  flow  region  of  the  Titan  IV  nozzle, 
because  of  overexpansion  of  the  flow  downstream  of  the  throat,  which  was  also 
observed  in  the  cold  flow  test. 
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8.  ONE-PHASE  HOT  FLOW,  y  =  1.19 


Based  on  the  assumption  that  the  measured  cold  flow  Mach  number  at  the 
motor  inlet  plane  is  applicable  to  a  full-scale  hot  firing,  the  one-phase 
numerical  calculation  for  a  combustion  gas  (y  =  1.19)  at  chamber  temperature 
5900°R  and  chamber  pressure  570  psia  (78.5  sec  into  the  burn)  in  the  Titan  IV 
SRM  is  also  carried  out.  The  assumption  is  valid  from  the  following 
observations : 

a.  The  flow  is  insensitive  to  a  Reynolds  number  change  from  a  chamber 
pressure  of  45  psia  to  150  psia  in  the  cold  flow  test. 

b.  At  the  motor  inlet  plane,  the  one-dimensional  Mach  number  difference 
between  y  =  1.4  and  y  =  1.19  is  less  than  2.4%  at  the  same  area 
ratio  =  8 . 2C . 

c.  The  Mach  number  at  the  motor  inlet  plane  is  very  small;  the  axisym- 
metric  calculations  from  a  characteristic  formulation  for  the  inlet 
flow  in  Reference  7  indicates  that  the  difference  in  Mach  number  at 
the  motor  inlet  plane  for  gases  with  y  =  1.4  and  y  =  1.19  is 
negligible . 

The  converged  Titan  IV  one-phase  hot  flow  transonic  solution  tar.es  2168 
integration  steps  which  require  5  min,  7  sec  computation  time  on  a  CRAY  X-MP14 
machine.  The  computed  three-dimensional  nozzle  hot  flow  discharge  coefficient 
is  0.9850,  based  on  a  converged  three-dimensional  gas  mass  flow  rate  of 
3920.07  lb/sec  and  a  one-dimensional  isentropic  gas  mass  flow  rate  of  3979.64 
lb/sec  for  the  full-scale  Titan  IV  motor  with  eroded  throat  diameter  3.3958  ft 
at  78.5  sec  into  the  burn.  The  computed  Mach  numbers  at  the  throat  are  1.305, 
1.292,  1.266,  and  0.824  at  TDC ,  90°,  BDC,  and  centerline,  respectively,  for 
the  hot  flow  with  y  =  1.19. 

An  isometric  projection  of  a  three-dimensional  Mach  number  contour  is 
given  in  Figure  14.  In  general,  the  overall  three-dimensional  flow  structure 
for  the  hot  firing  with  y  =  1.19  is  similar  to  that  of  the  cold  flow  with 
y  =  1.4.  Changing  y  from  1.4  to  1.19  has  the  effect  of  decreasing  the 
Mach  number  at  supersonic  speed  and  increasing  the  Mach  number  at  subsonic 
speed.  Note  that  the  Titan  IV  SRM  nozzle  has  a  larger  and  flatter  throat, 
hence  a  lower  throat  wall  Mach  number,  than  that  of  the  Titan  III  SRM  nozzle 
discussed  in  Reference  7. 
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9.  TWO- PHASE  HOT  FLOW,  GAS  y  =  1.19 


The  two-phase  flow  property  data  for  the  Titan  IV  are  shown  in  Table  1. 

The  computed  results  from  the  one-phase  y  =  1.19  hot  firing  are  taken  to  be 
the  initial  guess  for  the  fully  coupled,  gas-particle,  two-phase  flow  analysis. 
For  the  two-phase  hot  flow  with  a  single  particle  size  of  6  p  radius  in  the 
flow  field,  the  converged  solution  with  a  very  small  convergence  criterion  of 
0.001%  requires  5869  integration  steps  and  40  min  computation  time  on  a  CRAV 
X-MP14  super  computer.  Increasing  the  convergence  criterion  to  0.005%  will 
significantly  reduce  the  required  computation  time  to  less  than  10  min.  The 
discharge  coefficient  is  found  to  be  1.1764,  based  on  a  computed  three- 
dimensional  total  mass  flow  rate  4681.59  lb/sec  (gas-phase  3348.29  lb/sec  + 
particle  phase  1  333.30  Ib/sec)  and  a  cne-d iniens ional  isentropio  gas  mass  flow¬ 
rate  3979.64  lb/sec.  The  discharge  coefficient  is  greater  than  one,  when  a 
two-phase  mass  flow  rate  is  compared  with  a  one-dimensional  gas  mass  flow¬ 
rate.  The  computed  gas-phase  Mach  numbers  at  the  throat  are  1.152,  1.096, 
1.092,  and  0.721  at  TDC ,  90°,  EEC,  and  centerline,  respectively,  for  the 
two-phase  hot  flow.  The  gas-phase  Mach  number  contours  for  the  two-phase  hot 
flow  are  given  in  Figures  15  and  16,  where  the  sonic  surfaces  for  both  one- 
and  two-phase  flows  are  shown  for  comparison.  With  the  presence  of  particles 
in  the  flow  region,  gas-phase  expansion  is  slowed  down  significantly.  The 
Mach  number  in  the  present  two-phase  study  is  evaluated  from  the  local  sonic 
speed  for  the  gas-phase  with  y  =  1.19.  The  Mach  number  will  be  approximately 
2%  higher  than  that  shown  in  this  report,  if  it  is  computed  from  a  local  sonic 
speed  corresponding  to  an  equilibrium  gas-particle  mixture  (EGPM)  of  y  =  1.152 
based  on  the  equation  shown  in  Reference  7. 


Table  1.  Titan  IV  Flow  Property  Data 


Gas  Phase 

C  =  2.68  KJ/Kg-°K(0.64  Btu/lb  -  °R) 

p  m 

u  .  =  8.88  x  10_5Pa-s(5.97  x  10~5 lb  /ft-s) 

rtl  m 

Pr  =  0.45 
a  =  0.664 
y  =  1.19 
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Partic le  Phase 

C.  =  1.38  KJ/Kg-°K(0. 33  Btu/lb  -  °R) 
j  m 

m.  =  3203.69  Kg/m3(200  lb  /ft3) 
j  m 

4>  =  28.8% 


Figures  17  through  20  summarize  and  compare  the  Mach  number  and  pressure 
distribution  on  the  wall  and  centerlines  of  motor  and  nozzle  flow  region  for 
one-phase  cold  flow,  one-phase  hot  flow,  and  two-phase  hot  flow.  Although  the 
difference  in  flow  variables  between  the  one-phase  cold  flow  (y  =  1.4)  and 
the  hot  flow  (y  =  1.19)  is  negligible  in  the  low  subsonic  motor  flow  region, 
it  cannot  be  ignored  in  the  transonic  and  supersonic  nozzle  flow  region.  The 
difference  is  more  pronounced  when  the  particles  are  introduced  into  the 
flow.  A  higher  pressure  load  on  the  nozzle  wall  from  the  hot  flow  than  that 
from  the  cold  flow  is  evident  from  the  results  of  analysis.  The  design  of 
solid  rocket  motor/nozzle  and  evaluation  of  motor  performance,  therefore, 
should  not  rely  entirely  on  the  data  obtained  from  the  cold  flow  test. 


MOTOR  XI  (feet) 


Figure  17.  Comparison  of  Mach  Number  Distributions  (motor) 
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NOZZLE  X2  (fee!)  NOZZLE  X2  (feet)  NOZZLE  X2  (fee!) 


Figure  18.  Comparison  of  Mach  Number  Distributions  (nozzle) 


MOTOR  XI  (feel) 


Figure  19,  Comparison  of  Pressure  Distributions  (motor) 


- 1-PHASE  (7  =  1.40) 

- 1-PHASE  (7  =  1.19) 

- 2-PHASE  (7  =  1.19) 


NOZZLE  X2  (feet)  NOZZLE  X2  (feet)  NOZZLE  X2  (feet) 

Figure  20.  Comparison  of  Pressure  Distributions  (nozzle) 

For  a  giant  motor  like  Titan  IV  (10  ft  dia,  113  ft  long),  it  is  important 
to  have  accurate  information  on  total  pressure  drop  from  the  forward  end  of 
the  motor  to  the  aft  closure.  In  a  motor  firing  test,  the  forward  end 
pressure  can  be  conveniently  measured  (usually  through  the  igniter  boss).  But 
aft  closure  total  pressure  is  usually  unknown,  even  though  it  is  the  aft 
closure  total  pressure  that  directly  affects  performance  and  payload  delivery 
capability  of  a  motor,  especially  during  ground  level  vehicle  liftoff.  An 
empirical  approach  has  been  taken  in  the  past  for  calculating  the  total 
pressure  loss  factor  with  only  limited  success.  A  cold  flow  test,  although 
providing  useful  information,  is  often  carried  out  in  a  benign  gas-only 
one-phase  flow  environment  and  is  limited  to  the  number  of  propellant  burn-back 
surfaces  that  can  be  simulated.  In  this  regard,  if  static  pressure  in  the 
Liquid  Injection  Thrust  Vector  Control  (LITVC)  ports,  located  in  the  exit  cone 
downstream  of  the  nozzle  throat,  is  monitored  during  full-scale  Titan  IV 
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static  tests,  then  numerical  calculation  based  on  the  present  analysis  can  be 
iteratively  applied  for  defining  the  aft  closure  pressure  field  corresponding 
to  the  measured  static  pressure  at  downstream  exit  cone  locations. 

The  gas-phase  velocity  vector  plot  on  the  symmetry  plane  is  presented  in 
Figure  21.  A  three-dimensional  isometric  projection  of  the  gas-phase  velocity 
vector  plot  is  given  in  Figure  22.  These  velocity  plots  illustrate  that  the 
gas  flow  aligns  itself  from  the  chamber  toward  the  nozzle  and  then  follows  the 
general  trend  of  canting  direction. 


Figure  21.  Gas  Velocity  Vectors  on  Symmetry  Plane  (two-phase,  y  =  1.19) 


Figure  22.  Three-Dimensional  Gas  Velocity  Vectors  (two-phase,  y  =  1.19) 

The  three-dimensional  feature  of  two-phase  flow  in  the  Titan  IV  SRM  is 
portrayed  in  the  particle  density  contours  and  particle  velocity  vector  plots 
of  Figures  23  through  26.  Since  a  particle  flow  path  cannot  faithfully  follow 
a  steep  change  in  wall  slope,  particles  are  clustered  near  the  wall  with  a 
steep  entrance  angle.  In  the  nozzle  flow  region,  particle  density  on  the 
wall  differs  greatly  from  TDC  to  BDC  plane,  revealing  the  influence  of  non- 
axisymmetric  boundary  geometry  on  particle  flow  path.  A  distinctive  particle- 
free  zone  appears  in  the  calculated  results.  On  the  TDC  wall  with  a  steep 
slope  change,  the  particle-free  zone  occurs  upstream  of  the  throat.  On  the 
BDC  wall  with  a  gradual  slope  change,  the  particle-free  zone  occurs  down¬ 
stream  of  the  throat.  Moreover,  since  heavy  particles  cannot  effectively  turn 
around  the  corner  of  the  nozzle  entrance  wall  with  a  steep  slope  change  and 
tend  to  cluster  near  the  centerline  of  the  nozzle,  significant  reduction  in 
gas-phase  velocity  in  the  two-phase  flow  from  that  in  the  one-phase  flow  at 
the  centerline  near  the  nozzle  exit  plane  is  observed. 
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Figure  25.  Particle  Velocity  Vectors  on  Symmetry  Plane 
(two-phase,  y  =  1.19) 


FREE  ZONE 


Figure  26.  Three-Dimensional  Particle  Velocity  Vectors 
(two-phase,  y  =  1.19) 
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The  particle  velocity  lag  X^  and  gas-to-particle  temperature  ratio  X^, 
on  the  nozzle  throat  plane  and  the  nozzle  exit  plane  are  given  in  Figure  27. 

On  the  nozzle  throat  plane,  X^  and  X^,  are  higher  at  the  centerline  than 
those  near  the  wall.  On  the  nozzle  exit  plane,  the  reverse  is  true  because  of 
gas-phase  recompression  near  the  wall  downstream  of  the  throat,  which  reduces 
the  magnitude  of  relative  speed  between  gas-phase  and  particle-phase.  Cluster¬ 
ing  of  particles  near  the  nozzle  centerline  downstream  of  the  throat  also 
causes  reduction  in  gas-phase  speed,  resulting  in  a  slight  increase  of  X  from 
nozzle  throat  to  exit  plane.  Incidentally,  the  particle  Reynolds  number, 
which  is  a  measure  of  the  magnitude  of  relative  speed  between  gas-phase  and 
particle-phase,  varies  from  0  on  the  motor  entrance  plane  to  17  (TDC),  20  (90°), 
21  (BDC),  15  (centerline)  on  the  nozzle  entrance  plane,  to  295  (TDC),  31  (90°), 

31  (BDC),  20  (centerline)  on  the  nozzle  throat  plane  and  to  2 47  (TDC),  233  (90°), 
235  (BDC),  18  (centerline)  on  the  nozzle  exit  plane.  In  the  particle-free 
zone,  the  particle  velocity  is  zero  and  the  particle  Reynolds  number  becomes 
the  Reynolds  number  based  on  the  gas  speed  and  a  reference  particle  radius. 


NOZZLE  R2  (feet) 


Figure  27.  Particle  Velocity  Lag  and  Temperature  Ratio  (nozzle) 
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The  particle-phase  flow  fields  presented  in  Figures  23  through  27  are 
useful  to  rocket  nozzle  designers  for  particle  impact  erosion  calculations. 
Insulation  materials  such  as  silica  fiber  and  asbestos  filler-loaded  nitrile 
butadiene  rubber  (NBR)  aft  closure  insulation,  graphite-phenolic  nozzle 
throat,  and  silica-phenolic  exit  cone  and  exit  cone  extension  used  in  the 
Titan  IV  SRM  provide  thermal  protection  for  motor  metallic  structural 
components.  Adequate  insulation  design  without  undue  weight  penalty  is 
important  to  vehicle  survivability  in  a  hostile  launching  environment.  An 
investigation  of  enhanced  insulation  erosion  from  particle  impact  usually 
requires  knowledge  of  incident  particle  density,  flow  angle,  and  velocity. 

This  information  can  be  obtained  directly  from  a  fully  coupled,  two-phase 
analysis  such  as  the  one  shown  in  this  report. 

Finally,  a  post-processor  was  recently  developed  for  automated  color 
display  of  the  computational  results  from  the  present  analysis.  Figures  28 
and  29  compare  the  color  three-dimensional  Mach  number  contours  for  the 
gas-only  one-phase  and  fully  coupled  gas-particle  two-phase  flows  inside  the 
Titan  IV  SRM.  The  lines  in  these  figures  are  portions  of  the  interior 
computational  grid.  The  numbers  next  to  the  color  chart  indicate  contour 
values.  As  mentioned  previously,  the  Mach  number  for  a  two-phase  flow  is 
lower  than  for  a  one-phase  flow  and  can  be  seen  from  these  figures.  The  color 
three-dimensional  presssure  contours  (Fig.  30)  show  that  the  pressure  remains 
essentially  unchanged  on  most  of  the  aft  closure  interior  surface  and  drops 
rapidly  in  the  throat  region  of  the  nozzle.  The  overall  pressure  load  on  the 
aft  closure  and  nozzle  is  higher  for  the  two-phase  flow  than  for  the  one-phase 
flow.  The  pressure  distribution  shown  here  is  required  for  motor  performance 
and  thermostructural  analysis.  The  color  three-dimensional  particle  density 
contours  (Fig.  31)  reveal  high  concentration  of  particles  near  concave  regions 
of  the  aft  closure,  indicating  high  local  insulation  erosion.  This  has  been 
confirmed  by  the  data  from  two  full-scale  Titan  IV  static  motor  tests  and 
partially  explains  the  burn-through  failure  in  this  area  of  the  first 
seven-segment  MOL  test  motor  19  years  ago. 
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Figure  30.  Three-Dimensional  Gas  Pressure  Contour,  Titan  IV  Solid  Rocket 
Motor  (two-phase  flow,  v  =  1.19) 
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Figure  31.  Three-Dimensional  Particle  Density  Contour,  Titan  IV  Solid  Rocket 
Motor  (two-phase  flow,  y  =  1.19) 


10.  CONCLUSION’S 


In  this  study,  the  treatment  of  complicated  three-dimensional  motor  and 
nozzle  geometries  is  facilitated  through  use  of  three-dimensional,  body-fitted 
grid  systems.  Difficulty  and  confusion  associated  with  a  canted  nozzle  due  to 
different  sets  of  reference  coordinates  are  eliminated  through  introduction  of 
multiple  flow  regions  which  allow  precise  interpretation  of  flow  variables  in 
each  region.  The  idealized  fully  coupled,  gas-particle,  two-phase  flow 
governing  equations  in  three-dimensional  space  are  solved  numerically.  Cold 
flow  test  data  are  utilized  for  verification  of  the  operational  computer 
program  which  is  developed  from  a  careful  numerical  consideration.  Results  of 
the  analysis  are  presented  for  a  rocket  motor  of  practical  importance  in  space 
exploration.  Important  features  of  compressible  nozzle  flows  are  discussed; 
namely,  gas-only,  one-phase  flows  with  different  ratios  of  specific  heats  and 
a  fully  coupled,  gas-particle,  two-phase  flow  inside  a  canted  three-dimensional 
Titan  IV  SRM.  The  groundwork  laid  herein  will  be  helpful  to  the  further  study 
of  viscous,  three-dimensional,  multi-particle  sized,  two-phase  nozzle  flows. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer"  for 
national  security  projects,  specializing  in  advanced  military  space  systems. 
Providing  research  support,  the  corporation's  laboratory  Operations  conducts 
experimental  and  theoretical  investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  is  the  technical  staff's  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 

Aerophyslcs  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contami nation,  thermal  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  excimer  laser  development  including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out-of -f ield-of -view  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications; 
microwave  semiconductor  devices,  raicrowave/millimeter  wave  measurements, 
diagnostics  and  radiometry,  microwave/millimeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy-induced 
environments. 

Space  Sciences  Laboratory:  Magnetospheric ,  auroral  and  cosmic  ray 
physics,  wave-particle  Interactions,  magnetospheric  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomy, 
infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  syscems;  space 
instrumentation. 


